NPS ARCHIVE 
1969 

SMYTH, N. 



THE CROSS SECTION FOR THE FORMATION OF 
IN THE REACTION OF FAST PROTONS 
WITH METHANE 



I 



by 

Norman Robert Anton Smyth 



+ CVJ 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 



United States 

Naval Postgraduate School 




THESIS 



THE CROSS SECTION FOR THE FORMATION OF 
IN THE REACTION OF FAST PROTONS WITH METHANE 



by 



Norman Robert Anton Smyth 



June 1969 



T/u4 documj/U koA been approved {^OK pubtic Ke~ 
teoAe and AaJbi; -cU dut/iibution unloncted. 




The Cross Section for the Formation of 
in the Reaction of Fast Protons with Methane 



by 



Norman Robert Anton Smyth 
Captain, Canadian Armed Forces 
B.Sc., University of Alberta, I963 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1969 







ABSTRACT 

. + 

The capture cross section for the formation of m the 

+ + , . 
reaction H + CH# ”* + CH^ was measured at incident proton 

energies of 70, 85, 100, 150 and 200 eV and covering the scat- 
tering angles of 43 ° to 49 - 5 ° (lab coordinates). At 100 eV and 
below the curve of the cross section versus angle shows a sharp 
peak at about 46° whose position approaches the theoretical limit 
of 46.9° with increasing energy. Above 100 eV the peak was too 
small to be observed and only an upper limit can be placed on 

the value of the cross section. Typical values of the total 

-21 2 -22 
cross section are 2.0 x 10 cm at 70 eV and 7*6 x 10 at 

100 eV . The magnitude and energy dependence of the cross section 

as well as the angular position of the peak all are in essential 

agreement with the classical theory of ion-molecule rearrangement 

collisions proposed by Bates, Cook and Smith. 
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I. 



INTRODUCTION 



The correct description of complex atomic and molecular inter- 
actions is one of the major achievements of modern quantum theory. 
While there appears to be little doubt about the correctness of 
the theory in principle, it has become evident in practice that 
the mathematical complexity is such that only the very simplest 
cases can be treated with anything approaching complete rigor. 
Therefore many approximations are employed to simplify these 
mathematical computations. The question then, of the accuracy, 
nature and range of applicability of these approximations is of 
considerable practical importance. In general these questions can 
only be answered by comparison of the experimental results and the 
pr edictions of a specific approximate calculation. 

This paper presents the results of a series of measurements 
which confirm the predictions of a classical theory of ion- 
molecule rearrangement collisions at high impact energies as pro- 
posed by Bates, Cook and Smith [ l] . This theory uses an impulse 
type approximation to obtain the cross section for the capture of 
a light atom or ion from a target molecule by a fast moving pro- 
jectile. Applying this theory to the reaction: 

+ CH^ ^2 1 ) 

the cross section is predicted to show a sharp peak at 46.9^ with 

an upper limit on the magnitude of the cross section of 
/ -2 0 2 

1.4 X 10 cm at 100 eV incident proton energy. It is further 
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predicted that the magnitude of the cross section should decrease 

rapidly with increasing proton energy, approaching an energy de- 

-5.5 

pendence of E ’ asymptomatically at high energies (above 500 eV) . 
The range of validity of the theory extends from about 50 to 800 eV 
when applied to reaction (1). 

Our measurements were carried out at 70, 85? 100, 150, and 
200 eV, covering the angles of scatter from 43^ to 49.5^. We 
found essential agreement between the predictions of the theory 
and the experimental data, indicating that the classical impulse 
approximation proposed by Bates, Cook and Smith is a valid model 
for reaction (1) in the energy region considered. 
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II. THEORY 



A. CLASSICAL THEORY OF ION-MOLECULE REARRANGEMENT COLLISIONS 
A classical theory of ion-molecule rearrangement collisions 
at high impact energies has been proposed by Bates, Cook and 
Smith [ ll . The impact energies are assumed to be high enough 
so polarization forces and chemical binding energies of the 
colliding molecules can be ignored. The theory applies to various 
ion molecule rearrangement reactions of the type: 

1 k — 1 ' k ^ ' ' 

where X^ and Y are simple atoms or ions. The bar indicates which 
ions are fast in the laboratory coordinate system. An example of 
a rearrangement reaction of this type is: 

H*^ + CH^ H^ + CH^ ( 2 ) 

The basic assumption made by Bates, Cook and Smith, is that 
process (1) may be described by a classical impulse approximation 
similar to that developed by Thomas [2], for the description of 
electron capture, where each composite system is regarded as a 
loose cluster of atoms and ions. ( Eg H^ = H + h"*" and 
CH^ = C + 3H + H) . 

According to this model, we can consider the reaction to occur 
in the following manner: consider a particle of mass M^ moving with 
velocity v^ through molecules of loosely bound particles of mass 
M^ and M^^. Figure 1 shows the sequence of collisions that lead to 
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FIGURE I THE CAPTURE MECHANISM 
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the capture of by . First mass collides with mass • 

—4 

This causes to be scattered at an angle Qj and velocity v| 

— ♦ 

while recoils at an angle 0^ with velocity v^ . Mass then 
suffers a second binary collision with mass This causes 

—4 

to scatter at an angle 0^ with velocity v^^. 



Now if: (a) 






^2 



or vJ = v" 



And if: (b) 



^2 



(3) 



then masses and may have a relative energy of motion below 
that required for separation, and so, the two particles may com- 



bine and proceed on as one. If the particle with mass is a 



composite particle then its disruption can be avoided if 



(A) 



In the particular case of: 

«2 = “l »k 

condition (3) is satisfied if 0j ^ 45^ s^nd 0” ^ 90^ (see Section 

II of this Chapter), but then condition ( 4 ) is violated. In the 

case of the reaction: 

+ CH, + CH_ 

4 3 

where M = M , = M and M = M , the violation of condition 

1 H+ H k ^^3 

(4) is not relevant since is not a composite particle. 

Consider now, the first binary collision between and s^s 
shown in Fig. 2. The probability of scattering into the solid 
angle dn(0^)with velocity between v^ and v^ + dv^ is: 
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FIGURE 2. THE FIRST BINARY COLLISION 




FIGURE 3. THE SECOND BINARY COLLISION 
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where appropriate differential scattering cross 

section. Since the solid angle is the cone of semi-^angle 0^ within 
d0^ then: 

dQ(O^) = 2rr sin 0^ d0^ 
and equation (5) becomes: 

q = sin 0^ d0^ . (6) 

Consider now, the second binary collision between s^nd as 
shown in Fig. 3- The probability of having an impact parameter 
between ^ and ^ + d(^ and an azimuthal angle between ij) and 0 + d0 
at a distance r from is given by: 

p = ede — ^ (7) 

4tt r 

where r is the ^2”^k internuclear separation. But 

particles incident on with impact parameter p and azimuthal 
angle 0 are scattered into the solid angle dr(0”) at 0^. 

Hence: ^ ^ ®2 ^ dQ(0^) 



where appropriate differential scattering cross 

section . 

But: dO(0”) = sin 0^ d0” 



Hence equation (7) becomes; 



= ^ 2 k (® 2 ’) ®2 ^®2 



d0 

, ; 

4tt r 



( 8 ) 



Now, for to capture 
and must be contained within a 



- v” the relative velocity of 
volume of velocity space 
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determined by D, the mutual affinity of these two systems » In 
the high velocity limit. Bates, Cook and Smith [ l] , state this 
requirement as: 

3/2 



(v")^ sin 0" dO" d»/j dv" = [~ 



(9) 



where \ ^ and u = ^^2 case of 

Substituting equation (9) into equation (8) we get 

3/2 , 



P = 



2k 



3r 



2 ^ u, ^ 



(vp dv^ 



( 10 ) 



( 11 ) 



Now, the capture cross section arising from the 



sequence of binary collisions is: 



Q(M^-M 2 -Mj^) = Vpq 



( 12 ) 



where V is a dimensionless factor, less than, or equal to unity, 
which allows for the fact that and niay approach each other 
such that in the resulting state, the affinity is less than that 
assumed. 

Substituting equations (6) and (11) into (12) one gets: 

2tt cr^ ^(©’ )cr^, (0”) 

12^ 2^ 2k^ 2^ 



q(Mj^-M2"Mj^) = y\ 



3r 



/2D \3/2 

®2 (^3) 



Changing from Lab, to Center of Mass coordinates, one gets for 

the case of M, = « M, 

12 k 



Q(m^-M2~Mj^) = y- 



16tt a^2(90°) 

3r^ 




( 14 ) 



where and denote the differential scattering cross 



sections in the center of mass coordinate systems. This 
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differential cross section is at an energy of relative motion of 
2 

\ ^2^1 affinity D is at an M2"^k distance 

of r. To obtain the total cross section for capture, Q, one 

now sums over all binary collision sequences which lead to the 
same final capture process. For example, the total cross section 
for process (2) is 

Q = 4Q(H-H-C) 4Q(H-H*^-C)- 

This represents the two possible collision sequences shown below: 




Solution of Rutherford’s expression for the differential 
scattering cross section when v^ is high enough for the atomic 
field to be Couloumbic, shows that: 

Q(M^-M2-Mj^) ~ . (15) 

The range of validity of (l4) does not extend indefinitely 
as v^ increases. Two requirements must be satisfied: 

(i) The De* Broglie wavelength limitations, 

(ii) The Heisenberg uncertainty relations. 
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Bates, Cook and Smith [l], computed the required differential 
scattering cross sections using classical 

methods and the analytic representation of the relevant Hartree 
potentials as given by Byatt [ 3l • From the results of these 
computations listed in Ref. 1, the cross section for the re- 
arrangement reaction: 

+ CH, ^ + CH^ 

~ 4—23 

is determined to be as shown in Fig. 4* 



B. THE + CH, H^ + CH^ REACTION KINETICS 

~ 4—2 3 

The rearrangement mechanism for the reaction is shown in 

. + 

Fig. 1, where represents the H , the H and represents 

the ^^2 * determine the exact value of the scattering angle 

one has to consider the kinematics of the arrangement process. 

If the speed v^ is assumed to be large enough so the binding 

energy of M to M, can be ignored, then the conservation laws 

^ K 

can be applied in the usual manner . 

Consider the first binary collision of the reaction shown in 
Fig. 1. From conservation of energy and momentum we have: 






(16) 



M V 
1 1 



cos 0^ + cos 0^ 



( 17 ) 



0 = 



MiV| sin Q* - ^ 2^2 



(18) 



Now : 



M = M = M 
1 2 H 



Hence: the above equations become 



2 ,1x2 / , x2 

= (v{) + (vp 



( 19 ) 
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FIGURE 4 

PREDICTED CROSS SECTION FOR 

H * ♦ CH. — ► H ♦ ♦ CH, 

4 2 3 

FROM BATES COOK SMITH' 
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vJ cos 0' + v' cos 0' 
1 12 2 



( 20 ) 



V, = 



0 = v' 



s in 






s in 



®2 



( 21 ) 



Now 0^ can be eliminated between (20) and (21) by rearranging, 
squaring, and adding the two equations to yield: 

+ (v|)^- 2v^v| cos 0j = (v^)^ . 

But from (19): ~ '^l ” 



Thus : 


2 


(v|)^- 


2v,v* cos 




1 


1 1 


1 


Hence : 




(V{)^ 


= V^v| cos 





v| = cos 0| (22) 



Substituting equation (22) back into equation (19) we get: 

2 2 , ,2 , . ,2 
Vi = (cos 0|) + (v^) 

or: (1-cos^ 0^) = (v^)^ 

= v^sin 0^ ( 23 ) 

Equations (22) and (23) give us v| and in terms of and 
0| . Now we wish to obtain and in terms of and 0^* Thus 
we arrange equation (20) and (21) such that when we square and 
add them, we eliminate 0| and so obtain: 

^1 ^^ 2 ^^ ” ^^ 1^2 ^2 ” (^ i )^ • 



22 



But from (19) we have: 



2 2 2 
(Vi) = - (vp 

2 2 2 2 
Thus: + (v^) - cos 0^ = - (v^) 

Hence: cos 0^ . (24) 

Substituting equation (24) back into equation (19) gives: 

2 / , x2 2 2 

Vi = (v|) + cos 0^ 

or: v^(l-cos^0^) = (v|)^ 

v| = sin 0^ . ( 25 ) 

Consider now, the second binary collision. Again ignoring the 
binding energies we can write the conservation of energy and 
momentum equations as: 

% M2(v^)^ = % M2(v^')^ + ^ Mr ( 26 ) 





^2'’2 ~ ^2'^2" ^2 ^ ®k 


(27) 




0 = M„v" sin 0" - M, v, sin 0, 
2 2 2 k k k 


( 28 ) 


Now : 


Ml = M2 = 




and: 


^k = ^CH3 




Thus 


the above three equations reduce to: 






M 

/ 1x2 , „v2 ''k 2 

( V M = ( v*M + — V 

^ 2^ ^ 2^ k 

H 


(29) 




^k 

(V') = v^ cos 0M + _ cos 0J^ 

H 


(30) 
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0 = sin sin 0^^ 

H 



( 31 ) 



As we are not interested in 0, we can eliminate it by 

k 

suitably rearranging equations (3^) sind ( 31 ) 2(nd then squaring 
and adding. Thus we get: 



+ (v^)2 - 2v' v^-- cos Q'^=(~ vj , ( 32 ) 




But from equation (29) 



”k 2 _ 2 2 

r “k - ('’!> - (''2 > ■ 

H 



Thus (32) becomes 



(V^)^ + (v^)2 



2 '^k 2 

2v^v^ cos 0^ = — (vp - — (V-) 

H H 



or 



M, 



M, 



iT )■ ^ '^2^2 ®2 = ° 



H 



M, 



H 



M. 



Multiplying through by: 



H 



^k T2 



TTT yields: 



M -M / v" 
H k / 2 

M +M 
H k 




a 

V ' 



M 



H 



Mn-Mk 



COS 0^ 



= 0 



(33) 



M 



Now, let: a = 



H 






and b = 



^k-^H 



Then (33) becomes: 

2 




- 2 



2 / 



a cos 0^ 



b = 0 



( 34 ) 
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But now, for capture, we require that: 0^ -5:^ 0| + 0^ 



thus: cos 0” ^ cos (0| + 0^) 



ie. , cos 0” ^ cos 0^ cos 0^ - sin 0| sin 0^ amd using (2^) and 



( 25 ) to eliminate 0^ we have: 



vj 



V J 



. 2 \ 1 
cos 0" ^ cos 0* ( — I - sin 0’ I — 

2 ^ ^ V 



Hence equation ( 34 ) becomes: 




- 2a 



V ^ * V * 

2 3 

V * V 
2 1 



cos 0| 



v” V * 

2 1 . Of 

— r — Sin 0’ 
V ’ V 1 

2 1 



- b = 0 



or: (vp - 2a 



M,, I 



V”v 



2 2 



V”V 



cos 






2 1 



sin 0| 



- b(vp‘ 



= 0 . (35) 



But now, the second capture criteria is that: ^ 



But, from equation (22) we have: v^ = v^ cos 0^. 

Hence, we require for capture: v^ -5:; v^ cos 0^, (36) 

Substituting equations (22) and (36) into (35) we get: 

2 2 2 2 2 
V cos 0’ - 2a v’ [ V ’ cos 0*- v s in 0’ cos 0’1“b(v’) = 0. 

^ w ^ J. J. ^ ^ w 



Using (23) to eliminate v^ from this 

2 2 2 
v^cos 0| - 2av^sin0^ [v^sin0^cos 0|- 

2 2 2 
which reduces to: v^ cos 0| - b v^ 



equation we get: 

2 2 2 

v^sinO^cos 0|]“bv^sin 0^ 
sin^ = 0. 



0 



Hence , 



for 



capture to 
tan 0| = 



occur we must have: 




(37) 
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Now M, = M = 15*03506 and M = 1.00797 amu . Thus the angle at 
k. H 

3 

which capture occurs is: 

0| = 46.926° (38) 

Hence we expect the capture cross section to show a pronounced 

peak at this angle. From Fig. the magnitude of the cross 

; -20 2 

section IS predicted to be approximately 1.4 x 10 cm if the 
energy of the incident proton beam is 100 eV . 

C. THE RANGE OF VALIDITY OF THE CLASSICAL DESCRIPTION 

A lower limit on the range of validity of the classical 
treatment of the rearrangement collision is imposed by the 
assumption made in the theoretical development that the energies 
involved in the collision are much larger than the binding 
energies. Hence, when the energy of the incident proton beam E^ 
is less than, or equal to about 10 times the binding energy we 
expect the peak in the cross section to become less pronounced 
and more spread out. Thus the theory is valid when E^ ^ 50 eV . 

An upper limit on the range of validity of the theory is 
determined by the two requirements: 

(i) The Di Broghlie Wavelength A. of any particle involved 
must be much less than the smallest distance involved 
in the collision S. 

Hence we require: X « s (39) 

(ii) The uncertainty in energy introduced by specifying the 

transverse position of the proton must be much less than 
the binding energy of the final product. Now, from (39) 
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we know the uncertainty in the proton position X is less 
than the distance of closest approach S. 



But we can write the uncertainty relation: 



AxAp = 2 



as 



SAP » 5 



Thus 



^ 



or 



.2 

AE = - 

8m„S 

H 



2 • 



( 40 ) 



(41) 

(42) 



We then demand that the binding energy D be greater than AE, or; 



D 




(43) 



To find what restrictions equations (39) and (43) place on 
the energy of the incident proton beam we have to relate the 
distance of closest approach S to the energy of the proton beam 
. This can be done quite readily in the case of coulomb 
collisions. By considering the collisions to be coulomb it will 



be shown that: _ _1 



and 






fT 



where and are constants. Because S decreases faster than X 
as increases we see that equation (39) will indeed place an 
upper limit on the energy E^ . 
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In a Coulomb collision the point of closest approach is 



yielded by solving the equation: 



(see equation 3-4-3 of Ref. 4) 



(44) 



U, , the Center of Mass energy, can be expressed as U, = ^ ll v 
1 ^ 1 r el 

where u is the reduced mass and v , is the relative speed of 

r el ^ 

approach. 0{S), the potential energy at the distance of closest 
approach, is: 

0(S) = 



The impact parameter b is given by: 






(Ze)(Z’e') ^ 

b = ^ i- cot — 

^ ^rel 



(45) 



(see equation 3-8-8 of Ref. 4) 

Where (R) is the Center of Mass scattering angle. 

Now we have two collisions to consider: 

(i) In the first binary collision (H -H) we have (see Fig. 1) 

M 

® = 90 °, ^ (46) 

(ii) In the second binary collision (H-C) we have (see Fig. 1) 



(g) ^ QM = 90 °, u ^ 



el 



(47) 



and since ^ ^ E, we have vl ^ v^ . 

2 1 2 1 

Making these substitutions in (45) we get: 

2 



h'*'- H 



e^ 

E. 



"^H-C = 12 |- 



(48) 

(49) 
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wher e: 



E 



1 




2 



V 



1 * 



Hence for the H collision, equation (45) yields 




Which gives: S(h'^- H) 2.4 ~ . 



( 50 ) 



Similarily, for the H-C collision we get: 



12e 






which gives: 



s(H-c) ^ 7.-4 ^ . 

^1 



(51) 



Hence the H - H collision has the smallest distance of closest 
approach and hence equation ( 50 ) will establish the upper limits 
of validity of the classical approximation. 

Consider first the restrictions imposed by equation (39)- We 
require that: X « S 



or : 



— « 2.4 

1 

. E « 2.32 X 10 

• • JL 



Ergs 



or 



El « 



lc45 X 10 ^ eV. 



(52) 



, 4 

Thus the wavelength considerations require « 1-45 x 10 eV if 
our classical approximation is to be valid. 

Consider now the restrictions imposed by equation (43)* We 
require that: .2 

8Hh S 
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or u — - T • 

8m (2.4) 

H 

The binding energy of is 2.6 eV. Thus we get: 

El « ijSHii D (^4^) (53) 

= 3-93 X 10 ^ Ergs 

= 2.46 X 10*^^ eV. 



Hence equations ( 39) and ( 43) satisfied if the energy of 

the incident proton beam is less than 1,000 electron volts. Thus 
the classical description of process (2) is valid in the 75 to 
200 eV range where the experiment was performed. 
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III. EXPERIMENTAL APPARATUS 



The experimental apparatus used in this experiment is the 
result of the work of many peopleo The overall apparatus, as 
shown in Fig» 5y was assembled by Bush [5l> who demonstrated its 
satisfactory operation. 

A reasonably mono-energetic beam of hydrogen ions ( AE 2 eV ) 
is produced in the Duoplasma tr on , which was built and investigated 
by Carter [63<. The beam is mass analyzed in the mass spectrometer 
constructed by Strohshal [7]? sind ions proceed on to the scat- 
tering chamber which contains the target gas, The scattering 

cell is placed on the magnetic axis of a large cylindrical magnet 

at the position of maximum field. ions, produced by the col- 

+ 

lision of H with CH^, exit from the scattering cell and are 
focused by the axially symmetric non-uniform field of the focusing 
magnet to a point on the magnetic field axis, a distance from 
the scattering cello The detector which is located at Z , is an 

g 

electron multiplier capable of gains up to 10 . 

A . THE DUOPLASMATRON 

The duoplasmatr on is shown in Figo 6. With hydrogen in the 
source and filament power at 35 watts, one has about .1 amp of 
electron emission from the filament. An arc voltage of 400 volts 
between the filament and the Z*°electrode accelerates these elec- 
trons through the hydrogen gas, ionizing the hydrogen and forming 
a plasma. Once the plasma has been formed, an arc current of 
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FIGURE 5. SCHEMATIC OF THE EXPERIMENTAL APPARATUS. 
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1 ampere can be sustained between the filament and the Z-electrode 
by a voltage of about 250 volts. 

To extract the plasma, the positive arc voltage is switched 
from the 2-electrode to the anode, (see Fig. 7 ) which is itself at 
a voltage E above ground. The energy of the hydrogen ion beam thus 
extracted will be (Ee) for singly charged ions. 

To enhance the beam intensity a simple accel-decel system is 
employed. Figure 7 shows lens 1 at some arbitrary negative po- 
tential V. Thus the hydrogen ions are actually extracted from 
the duoplasma tr on at an energy (V+E)e. A ^’sock^* is fitted over 
the focusing and deflector lenses and attached to lens 1 so as to 
float at this negative potential V, and thus prevent the beam from 
^*seeing^* ground potential in this region. Thus the beam of hydro- 
gen ions passes through the focusing and deflector lenses at an 
energy (V+E)e. Once the beam leaves the region of the sock it 
^^sees^’ ground potential on the walls of the vacuum chamber and 
slows to an energy Ee and passes on to the mass analyzer . 

This simple accel-decel system increased the beam intensity by 

a factor of • Thus we were able to obtain an analyzed beam of 

100 eV protons (with diameter .5 cm at the scattering cell 80 cm 

-7 

from the source) which was slightly greater than 1 x 10 amps. 
Carter [6] has shown that the energy spread of the proton beam is 
less than 2.5%. 
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FIGURE 7. 

CIRCUIT DIAGRAM FOR DUOPLASMATRON 
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B. THE MASS ANALYZER 



The hydrogen ion beam is mass analyzed in the analyzer con- 
structed and calibrated by Strohsahl [7lo It employs a 45° 

bending angle (see Fig* 8) and has a resolution of — = 40 with 

Am 

beam transmission of about 90 %* 

Figure 9 shows the current through the mass analyzer magnet 
that will pass any particular species in the energy region from 
50 to 650 volts. Hence *45 amps allows a 100 eV beam of protons 

■f- + 

to travel into the scattering cell, whereas the H^ and H are not 
bent enough by the magnetic field to be transmitted (see Fig* 5) • 
The mass analyzer is carefully aligned with the scattering 
cell and the detector, all of which are positioned on the axis 
of the focusing magnet. (See section F below.) 



C. THE SCATTERING CELL 

The scattering cell, shown in Fig. 10, was designed and 

tested by Bush [l]. The target gas, methane, is bled into the 

scattering cell and is accurately controlled by a Variable Leak 

valve constructed by Granville Philips. The pressure of the CH^ 

-*3 

gas in the scattering chamber (approximately 10 torr ) is mea- 
sured by the VG lA ion gauge. The gas is of research grade, 

(pur ity 99 • 65 %) . 

The incident proton beam passes first through an adjustable 

aperature (generally set at .5 cm diameter) and then into the 

. + , + 
scattering chamber. Here the H reacts with the CH^ to form 

whose escape from the scattering cell is limited by the geometry 

of the exit slit. The exit slit is an annular opening of 3^0^, 
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SCALE INCHES 



Figure 8. Schematic of Mass ^ectrometer. 
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MAGNET CURRENT- AMPS 
FIGURE 9 

SPECTROMETER OPERATION CURVES 
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Figure 10. The Scattering Cell and Beam Collector 




reduced by support legs to about 290 ^, around the scattering cell, 
separating the front from the rear of the scattering cell. The 
exit wall of the front section of the scattering cell slopes up- 
wards at 49° and the exit wall of the rear section of the scat- 
tering cell slopes upward at 36^. Appendix III of Ref. 5 shows 
how these two angles determine the target thickness. By varying 
the separation of the front and rear sections of the scattering 
cell, the target thickness can be increased from 0 to 0.8 cm. 

The ions collected by the front and rear sections of the 
scattering cell are measured by separate ammeters. The beam col- 
lector collects the unscattered beam that passes through the 
scattering chamber. This current, plus the current from the rear 
section of the scattering cell, is recorded as the total proton 
current incident on the CH^ target molecules. 

Three grids on the beam collector are given potentials to 
suppress secondary electron emission and prevents slow background 
ions from reaching the collector plate. The beam collector and 
its grids are placed inside a metal cup which shields them from 
other charged particles and also acts as a gas cap over the rear 
of the scattering chamber. 

The beam collector assembly is mounted on a rod which passes to 
the outside of the vacuum system, allowing the assembly to be 
pulled away from the rear section of the scattering cell and swung 
out of the path of the incident ion beam. This allows alignment 
of the entire apparatus and also allows the incident beam to fall 
directly on the detector when desired for calibration purposes. 
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D. THE FOCUSING MAGNET 



The magnetic field of the f ecus ing magnet is nonuniform and 
axially symmetric and has the focusing property of bending the tra 
jectories of charged particles so they follow a path similar to 
that shown in Fig. 11. 

If one considers the simpler case of an axially symmetric 
uniform magnetic fields it can be shown that the trajectory of a 
charged particle is a helix (see Ref. 5)- The particle will cross 
the field line on which it originated after one complete revolu- 
tion. Then, as shown in Fig. 11, particles of charge q leaving 
the source with momentum p = mv at an angle Q will cross the 
magnetic axis at a distance given by: 



_ 2tt p cos Q 
qB 



( 1 ) 



where B is the magnetic field strength. Thus a detector located 
at will only sense these ions scattered out of the scattering 
cell placed at Z = 0 which have momentum p, charge q, and scat- 
tering angle ©. If one now holds the detector fixed at Z^ and 
increases B then particles of larger p cos 0 values will strike 
the detector. Alternately by holding B fixed and varying Z^ we 
can detect particles scattered from the source at different values 
of p cos 0, Since equation 1 is independent of azimuthal angle 0, 

the entire 3^0^ can be observed at Z . For most conventional scat 

o 

tering experiments A© ^ 1^. Hence, the solid angle is enhanced by 
a factor of approximately 3©© « 

In our experimental apparatus the magnetic field is NOT uni- 
form, as it is produced by a *^thin” solenoid shown in Fig. 5* 
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Figure 11. Orbit of Charged Particle in Uniform Magnetic Field 



This complicates the trajectory equations and makes numerical 

inter gr at ions necessary. The results of numerical inter gr at ions 
+ 

for the ion are shown in Fig. 12. However, equation 1 the 
result obtained for the uniform field case is qualitatively the 
same as in our NON~uniform field configuration. 

The magnetic field of the focusing magnet was measured by 
Kelly [8]. Figure 13 shows the axial component of the magnetic 
field on the magnetic axis. Figure l4 shows the values of the 
axial component of the magnetic field measured at various r and 
0 values at the center of the coil. Figure 15 ^ and b show the 
radial components of the magnetic field as measured with a Hall 
Probe(2^t = 34 cm) across the geometric axis of the coil, both 
vertically and horizontally. This demonstrates that the magnetic 
and geometric axis of the magnetic coil are aligned to within 
0.1%. This allows use of the coil geometric axis and the coil 
center -line as reference for alignment of the entire system along 
the magnetic field axis. 

E. THE SCATTERED ION DETECTOR 

The beam detector and amplifier is a Model 3^6 Magnetic Elec- 
tron Multiplier whose schematic is shown in Fig. l6 , The multi- 
plier is connected to a Bendix 1122 Power Supply as shown in 
Fig. 17* The output of the multiplier is fed to the input head 
of a Keithly 640 Vibrating Capaciter Electrometer. 

The Multiplier, developed by Goodrich and Wiley, provides re- 
producible current gains up to 10^ when operated at the potentials 
shown on Fig. 17* When operating as the detector of scattered 
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ions in our experimental apparatus, the multiplier has to be oper- 
ated in a magnetic field of up to 200 gauss. The effect of this 
magnetic field is to reduce the gain of the multiplier; increasing 
field, causes the gain to decrease. 

The multiplier is operated inside an aluminum shield to pro- 
tect it from random background. The various strip potentials were 
chosen to give the highest gain and the lowest noise. The multi- 
plier gain was measured at different magnetic field values in the 
manner listed in Appendix I and fitted by a 12th order polynomial. 
See Fig. l8 . The fitted polynomial was used in the actual cross 
section calculations . 

Bush, has shown that the gain does not vary with gas pressure 
-Zi 

up to 2,0 X 10 torr , a pressure much higher than the multiplier 
ever encounters when data is being taken. The gain is independant 
of the incident ion mass and of the incident current, provided 
the anode current is below its saturation value of 10 ^ amperes. 

To prevent slow ions from reaching the detector, a grid is 
placed over the detector shield aperature as shown in Fig. 19. 

The grid also permits a crude analysis of the scattered beam energy. 
See section A, Chapter IV. 

F. SYSTEM ALIGNMENT 

In order to detect particles scattered at a particular angle 
the detector must be placed on the focusing magnetic field axis at 
the point Z where the scattered particles focus on the axis. The 
program SOLANG which calculates the intercept distance assumes 
that the scattering took place on the axis at the center of the 
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Figure 19* Dector Aperature Geometry. 
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focusing magnetic field. Therefore accuracy in the angle of scat- 
ter determination depends on how well the scattering center and 
the detector are aligned to the magnetic field. 

Figure 15 shows that the geometric center line of the focusing 
magnetic coil is collinear with the magnetic field axis to within 
0.1%. This geometric centerline was used as the reference for 
aligning all components to the magnetic field axis. To perform 
the alignment, cross hairs were placed on either side of the mag- 
netic spool to mark its axis. Then a small laboratory laser was 
positioned so its beam lay on the magnetic axis established by the 
two crosshairs. The crosshairs were then removed and the line 
established by the laser beam was used as the reference for 
aligning all other components to the focusing magnet axis. This 
method made possible the alignment of the scattering cell, the 
detector, and the focusing lenses to the magnetic field axis. The 
mass analyzer was constructed by Strohsahl, (see Ref. 7 ) so that 
the ion beam, after analysis, travels down the geometric axis. 
Hence attaching crosshairs to the front and rear of the mass 
analyzer, it was adjusted into alignment with the laser beam. 
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IV . THE EQUATIONS DESCRIBING THE MOTION OF THE 



H^ ION AND THEIR SOLUTION 



A. THE TRAJECTORY EQUATION 

The equation of motion of the ion in the nonuniform axially 
symmetric magnetic field of the focusing magnet, (equation ( 5 ))? 
is expressed in terms of the magnetic vector potential A. There- 
fore the measured magnetic field B of the focusing magnet must be 
related to this quantity before the trajectory equation can be 
solved. The expansion of A in terms of the field B is shown by 
Bush [ 5l to be: 



Aj^(r ,z) 



= T. 



Izil 



n=o nl (n+ 1 ) I 



(|) 



2 n+l (r= 0 ,z) 

z ' ' 



bz 



2 n 



A^(r,z) = 0 
A^(r ,z) = 0 



( 1 ) 

( 2 ) 

( 3 ) 



Gagliano [ 9 ] fitted the actual magnetic field measurements along 
the r = 0 axis with a 12 th order polynomial in z such that: 

( 4 , 



where I is the current through the magnetic coil. Hence we know 
B (r = z) and A 0 (r,z) for any given focusing magnet current I. 
Gagliano, and subsequently. Bush, has shown that the equation of 
motion of a particle of charge q and mass m with momentum p in our 
field configuration is given by: 
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0 



( 5 ) 



d^r 



1 + 



fe) 



[k^ 



A^] - 



3? 



^0^ 



a A, 

^0 ar 



0 _ 



wher e 



‘^ = 4 

q 



( 6 ) 



To solve this second order nonlinear differential equation by 
means of the digital computer , it roust be reexpressed as two first 
order differential equations. This is done by rewriting equation 
(5) as follows: 



d 

dz 



dz 



dr 


'a N 


dz 


. ^0 az J 


, 2 
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■ ^0 



^1 /i'lt 

hz J ^ \dz) 1^0 Sz J 



2 2 
- A0 



a A, 



0 



"0 5r 

I 2 ^ 2 ' 

-A0 



Ktiz) 1^0 dr 



.2 ,2 



( 7 ) 



Define: 



DERY (1) 
DERY (2> 



= r 



dr 

dz 



a A, 



0 



FUNC ( 4 ) = 



FUNC ( 5 ) = 



”0 dz 

. 2 ^2 

- A0 



a A, 



0 



"0 ar 

, 2 «2 



Then equation ( 7 ) becomes 

[dERY(1)] = DERY(2) 



d 

dz 



( 8 ) 
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^ [dERY( 2)] = DERY(2) x FUNC(4) + [dERY(2)J^x FUNC(4) (9) 
- FUNC(5) - [dERY(2)J^x FUNC(5) 

We now have two first order nonlinear differential equations 
which can be solved by the use of the digital computer . 

From equation (6) we see we need to know the momentum of the 
scattered particle. But p = 2mE^ where is the energy of the 
scattered particle. Hence 



k 



2 




( 10 ) 



where is the energy of the scattered particle expressed in 
electron volts. Thus if E^ is known, equations (8) and (9) can 
be solved to yield r as a function of Z. 

B. THE ENERGY OF THE H^ 

^ + 

To determine E , the energy of the scattered H particle, 
s f ^ 

consider the following sketch of the kinematics of the reaction: 

H*^ + CH, ^ h!^ + CH^. 

— 4—2 3 



Before Collision After Collision 




where m represents the mass of the H , the incident ion, and the 
o ^ 

m^ represents the mass of the CH^, the target molecule, and m^ 

+ 

represents the mass of the the scattered ion, and the 
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4^ 
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^^0 



q 



= 0 (5) 



( 6 ) 



To solve this second order nonlinear differential equation by 
means of the digital computer , it must be reexpressed as two first 
order differential equations* This is done by rewriting equation 
(5) as follows: 
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Then equation (7) becomes 

[dERY(1)] = DERY(2) 
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dz 
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^ [dERY(2)] = DERY(2) x FUNC(4) + [dERY(2)]^x FUNC(4) (9) 
- FUNC(5) - [DERY(2)j^x FUNC(5) 



We now have two first order nonlinear differential equations 
which can be solved by the use of the digital computer , 

From equation (6) we see we need to know the momentum of the 
scattered particle. But p = ^”2m^ where is the energy of the 
scattered particle. Hence 



k 



2 




( 10 ) 



where is the energy of the scattered particle expressed in 
electron volts. Thus if E^ is known, equations (8) and (9) can 
be solved to yield r as a function of Z. 

B. THE ENERGY OF THE H^ 

To determine E , the. energy of, the scattered H particle, 
consider the following sketch of the kinematics of the reaction: 
h"^ + CH, ^ h!^ + CH^. 

- 4—23 



Before Collision After Collision 




where m^ represents the mass of the H , the incident ion, and the 

m^ represents the mass of the CH^, the target molecule, and m^ 

+ 

represents the mass of the H^ ? the scattered ion, and the m^ 
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represents the mass of the CH« , the remaining molecule. Ignoring 

j 

the binding energies, we can write the following: 

From the conservation of energy we have; 

2 2 2 

%mv =^mv +%mv (11) 

oo ss rr ^ ' 

From conservation of momentum we have: 



mV = m V cos 0 + m v cos 0 
o o s s r r 



( 12 ) 



0 = m V sin 0 - m v sin 0 
s s r r 



(13) 



We wish to solve these three equations to yield the relation; 

I". > I" » ® ) 

s s o r So 

Squaring (12) and (13) sind adding, we eliminate 0 and get: 

2 2 2 2 2 2 
m V - 2m V m V cos 0 + m v = m v 
oo ooss ss rr 

2 2 2 2 
But from (11) mv =mmv - mmv 
^ ^ rr roo rss 

Thus we have on substituting this into the above equation: 

2 2 2 2 2 2 
m V - 2m V m V cos 0+mv -mmv +mmv = 0 
oo ooss ss roo rss 



or 



2m E - 4 \1 m E m E 
o o y o o s 5 



cos 0 + 2m E ~2mE +2mE = 0 



s s 



r o r s 



Thus we get on rearranging: 

E(m +m)+E(ra - m)=2 \(m m E E ' cos 0 (l4) 

s' s 1 ' o o r ^ o s o s 



Now let: a = m^ + m^ and b = m^ - m^ 



( 15 ) 



Squaring (l4) we get: 
2 2 



2 2 r 2 1 2 2 

a E + 2ab E - 4 m m E cos 0 E + b E = 0 
sL o oso -'s o 
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Solving for E we obtain: 



E = 



2 2 2 2 2 2 
4 m m E cos 0 - 2abE + \ (2abE - 4 ni m E cos 0) - 4a b E 

oso o— V o oso ' o 



2a 



As the + root only yields physical results, in our case we have, 
on substituting equation (15) for a and b: 



E = E 
s o 



2mm COS0 + (m -m )(m + m ) 
o s ^ r o'^ ^ s r ^ 

(m + m ) 
s r ' 



2 2 4 



m m 



cos^0 - 4(ni )(m -m )m m cos 0 

•\ c V'' ■v‘ f 



( m + m ) ' 
s r ' 



Now: m = m ^ 

r + 



m — m . 
o 



m = m , — 2m 

S c 



Thus the above equation simplifies to: 



E = E 
s o 



2 2 2 2 
m + m m - 2m + 4ni cos 0 
r r o o o 



(m + 2m )' 
r o^ 



(16) 



^ 2 2 2 2 
2 m^ cos 0 \ m^ + m m - 2m _ + 2m cos 0 



r o 



( m + 2m ) 
' r o' 



Thus selecting 0 and E^ we can calculate E^ , the energy of the 

scattered H , from equation (l6). 

^ E 

s 

Figure 20, a plot of the ratio — versus the scattering angle 

E 

o 

0, shows that in the 10^ range from 4l^ to 51 j changes by less 
than 3%* 
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C. THE METHOD OF SOLUTION 



Assume protons of energy and a focusing magnetic current 

to have been selected. Then equation (4) yields the magnetic 

field B^(r = 0,z). Equation (l6) yields the energy of the H^ ^ls 

2 

a function of the scattering angle, 0. From equation (10) k is 

+ 

then known for the H^ and we are ready to solve equations (8) and 
(9). The solution is a trajectory similar to the ones shown in 
Fig. 7. The distance where r is zero is recorded. This is 
where the detector must be located to detect H^ particles scat- 
tered out of the scattering cell at the angle 0 when the focusing 
magnet current is set at and the energy of the incident beam 

of protons is E . 

^ o 

If the detector distance Z is changed and I and E are held 

o o o 

+ 

constant then H^ particles scattered out of the scattering cell 

at different scattering angles 0 are detected. 

Alternately, one can hold E and Z constant and vary I . This 

o o ' o 

also allows particles of different scattering angles to be detected. 
D. THE COMPUTER SOLUTION 

The computer solution of equations (8) and (9) was originally 
worked out by Gagliano, but has since been rewritten. Program 
SOLANG, explained in Appendix II, employs the Hamming Modified 
Pr edic tor -Cor r ector method to integrate equations (8) and (9). 

This is the DHPGC subroutine of the IBM System 36 O -67 Scientific 
Subr out in e Package. 
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The trajectory equations are solved for a range of scattering 
angles for a particular and The results are stored and 

later used to compute the solid angle subtended by the detector as 
it intercepts each of these trajectories. 

E. THE SOLID ANGLE 

The solid angle subtended by the detector is given by the 
relation: 

dO (Q,0) = sin 0 dO d0 

But since our system is axially symmetric this becomes 

dQ(0) = sin 0 dQj^d0 

2Q0 

^360^ ^ 

= 1.612 TT sin 0 dO. 

This is because only 290^ of the available 3^0^ of the axial angle 

0 are accessible to the detector due to the construction of the 

scattering cell. (See Chapter II, Section 3-) 

From the trajectory plots on Fig. 15 it can be seen that all 

particles scattered into 0 will cross the axis at Z within AZ . 

This AZ determines the angular resolution AQ which the detector 
o 

sees and so allows determination of the solid angle. But AZ is 

^ o 

in turn determined by the detector geometry and the particular 
trajectory under observation. The relation between AQ and AZ^ can 
be found from the functional relationship- 

= f(E, B, 0, , T], c) 

where E = energy of incident ion. 
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B = magnetic field 

^ = target width from system axis along a radial line 
Tl = target thickness from system center along the Z axis 
C = constant for each reaction including mass ratios, 
inelastic energy losses, etc. 

Differentiating this gives: 






dT ^ bT 



AB 



ae 



A0 



dZ dZ 

o o ^ 



5z 

o 

5C 



AC 



But AB = 0 and AC = 0 for any specific reaction and measurement. 

fdZ \ 



Thus, solving for AO we get: 



AZ - 
o 



AO = 






ATI 



E,Q, 



dZ 
( 

SO 









No analytical expression exists for the various differentials. 
However, computer trajectory results obtained by varying E,0,f^ 
and T1 show the effect of the differentials is to change AO by no 
more than by 10% (see Ref. 10). Thus the angular spread can be 
expressed as: 



AO 




AZ 

o 

E 



(17) 



Hence the solid angle becomes 



dQ (0) 



1 .6l2n sin 




(18) 



Now, since AZ^ depends on the angle 0^ (see Fig. 21 and 
Ref. 5 for further details) and 0^ in turn depends on the scat- 
tering angle 0 we see that dfi must be evaluated at each angle of 
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Figure 21. Angular and Energy Acceptance of Dector 
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scatter. This is implemented in the latter subroutine of the 
computer program SOLANG . 

Computer results indicate that the angular resolution given 
by equation (17) is approximately 1^. 

F. THE CROSS SECTION 

Consider a beam of monoener get ic ions incident on a target 
of molecules. A detector located at an angle Q subtends solid 
angle d 0 { Q ) . 




The number of particles scattered into the detector, , is 

usually measured as a current I where: 

s 

I = N dO (0) 
s s ' ' 



The cross 



s ection 



is then 
dJ (0) = 



def ined as : 
I 



s 

I . dfi(0) ' 
It ' 



(19) 



But I 

s 

the actual 



is magnified G times by 
detector current measur 



the multiplier 
ed is : 



ga in . 



Hence 
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Hence ( 19 ) becomes; 



I . G N d 0 ( 0 ) 

But the number of target particles, is given by; 

= P P t 

t SC 



Where = pressure of target gas in scattering cell in torr 



1 6 

^ = 3*536 X 10 par t icles/cm*^ torr 

t = target length in cm 



Hence we have: 



da( 0 ) = 



cm 



I. G P P t dfi(O) 

Because of the construction of the scattering cell (see 
Section 3 of Chapter II) the target thickness depends on the 
scattering angle 0 . Bush, has shown this relation to be: 

t = s - g ^ gO 

Sin 0 



( 20 ) 



t = S 



if 36° < 0 < 49° (21) 



t = s - ^ ^ 

s in 0 

where S is the minimum separation of the front of the scattering 
cell from the rear of the scattering cell, (which is adjustable 
from outside the vacuum system). 

Computer Program CSVSZ, listed in Appendix I, calculates the 
differential cross section in the Lab. and Center of Mass coordin- 
ate systems using equations (20) and (21) and the solid angle 
determined by the program SOLANG . 
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V. EXPERIMENTAL RESULTS 



A. ELIMINATION OF BACKGROUND 

When H collides with CH^ three general types of collision may 

occur which can produce large detector signals. 

Event (a) The rearrangement collision where the h"*" captures a 

+ 

proton to form • 

Event (b) A scattering collision of the h"*" with a hydrogen atom. 
Event (c) A scattering collision of the h"*" with a carbon atom. 
Events (b) and (c) are expected to have cross sections that are 
many orders of magnitude larger than those of event (a). Hence 
appropriate steps were taken to insure events (b) and (c) were 
not measured by our detector. To eliminate events (b) and (c) 
we look first at the energy and momentum of the scattered parti- 
cles in all three events. 

Consider first event (a). The capture process may be approxi- 
mated by the reaction shown below: 

+ CH, ^ + CH^ 

“ 4 “ 3 




The energy of the incident H particle is: 



E 

o 



1 2 
% m V 
o o 



2 




o 



( 1 ) 
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Now, from conservation of momentum we see the momentum of the 
+ 

scattered particle, represented by p^, related to p^ by; 



P 



o 



cos 45 




( 2 ) 



where we assume the mass of the represented by m^ , to be 

infinitely large. Hence in the capture process, the scattered 
particle has momentum p^ = >[2 p^. The energy of the scattered 
particle is then; 



E 

s 




+ . 

But for the particle, m = 2 m 
2 ^ s o 



Hence; 



E 

s 



2 




2m 

o 



E 

o 



( 3 ) 



Consider now event (b), the scattering of the h"*” from a 
hydrogen atom; 

h"^ + H h"^ + H 




Hence: 



Then: 



Pq ^ ^ Ps ^ 'P~ Pj 

_ Po 

2m 4m ^o- 

S O 



( 4 ) 

( 5 ) 
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Thus the elastically scattered from H has a momentum that is 

Po 

approximately — — and an energy that is approximately \ E . 

Finally consider event (c), the scattering of the from a 

carbon atom: ^ ^ 

H + C H + C 



V 




Assuming the carbon to be infinitely massive we see that; 

E = E (6) 

so ' ' 

2 2 

Hence: ^s _ ^o 

2m ” 2m 
s o 

But; m = m 

s o 

Thus: (7) 

Hence the elastically scattered h"*" from C has momentum p and 
energy E . 

Summarizing the results of this very approximate analysis in 

Table 1 we get for an incident beam of protons of energy E^ and 

momentum p the information shown. Note that the scattered H**^ 
o 

produced in events (b) and (c) has momentum less than the H^ 
produced in event (a). 



68 



Possible Events 



Event (a) Capture 
+ CH 

— 4—2 



Event (b) Scatter 
+ H + H 



Event (c) Scatter 
+ C h"^ + C 



Scattered 
Par t icle 



Momentum of 
Scattered Particle 



Energy of 
Scattered 
Particle 



3 



+ 




'F Po 



E 

o 




E 

o 



H 



TABLE 1. The Three Possible Scattering Processess. 

Plots of the scattered ion trajectories resulting from each 
of the three events tabulated are shown in Fig. 22, for the parti- 
cular case of E^ = 100 eV and I = 6.0 amps. Note that the two 

elastically scattered ions cross the axis at a Z less than that 

+ . + . 

of the ion. To prevent these unwanted H particles from 

reaching the detector, a baffle, 5 cm in radius was positioned 

over the beam collector as shown in Fig. 22. This allowed only 
+ 

the produced by the rearrangement collision to reach the 
detector . 

Figure 23 contains plots of the detector current versus the 
detector grid voltage for the case of E^ = 100 eV and I = 6.0 amps 
and with the detector located at Z = 46 cm (this corresponds to a 
scattering angle 0 of 46.92°). Figure 23 -(i) was recorded before 
the baffle was placed in the path of the ions and Fig. 23 - (ii) 
was recorded after the baffle was positioned, as shown in Fig. 22. 
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Consider first, Fig. 23 - (i). A grid voltage of 25 volts 
eliminates the slow ions. As the detector grid voltage is varied 
from 25 to 50 volts, we continue to collect the ions scattered from 
the scattering cell due to processess (a) and (b) and (c). Between 
60 and 90 volts we collect only these particles with energies 
greater than 90 volts; ie, the ions which result from processes 
(a) and (c) only. 

Figure 23-(ii) now shows that the baffle does indeed eliminate 
particles resulting from processes (b) and (c). Thus with the 
grid voltage set anywhere between 30 and 90 volts, and the baffel 
in position, we only detect the that results from the re- 
arrangement collision. 



B. DEMONSTRATION OF VALIDITY OF DETECTOR SIGNAL 



From equation (20) of Chapter IV we have: 
da(0) = :::: ^ , cm 



li ^ tG dQ{Q) 



( 8 ) 



If the angle of the scatter 0 is 
duces to: 



da 



K 



“D 



I .P 
1 sc 



held constant. 



then this re- 



(9) 



where K = constant. 

Consider now the case where I^ is held constant. Then (9) 
reduces to: 



D 1 sc 



( 10 ) 



Where C^ = constant. Experimental data plotted in Fig. 24 obeys 
this linear relationship. The fact that I^ does not extrapolate to 
zero is due to the background. 
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Dettctor current (Amps) 
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Consider next the case where P 



is held constant. Then 



sc 



equation (9) reduces to: 




( 11 ) 



where = constant. Experimental data plotted in Fig. 25 follows 
this linear relationship. Thus our detector signal has the proper 
dependance on the various parameters and behaves like a scattered 
particle current. 

For the differential scattering measurements made with 100 eV 

incident protons, the proton current at the beam collector was 
-7 

approximately 10 amps and the current to the rear of the scat- 

-8 

tering cell was about 10 amps. The methane gas pressure was 

-4 -3 

varied from 1 x 10 torr to 3 10 torr , The detector signal 

-12 

under these operating conditions varied from 1 x 10 amps to 
3 X 10 amps. This represents a flux of a few hundred particles 
per second. With no gas in the scattering chamber the detector 

-13 

current fell to about 5*0 x 10 amps. 

The UGIA Ion Gauge used to measure the methane pressure in the 

scattering cell was calibrated using a Capacitance Manometer. The 

relationship between the pressure measured by the ion gauge 

(P and the actual pressure which was measured by the mano 

measured^ 

meter (P , ) is linear in the region of interest as is shown 

^ actual^ 

below . 



2 . 0 -" 




slope = 0.625 



P 



actual X 10 
torr 



-3 



1.0 



1.0 



P 



measur ed 



2.0 



3.0 

X 10 "3 torr 
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-12 




Figure 24. Detector Current versus Pressure of Target Gas 



C3 




-8 

X 10 lamps , 



Figure 25. Detector Current versus Incident Current. 
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C. h'*’ + CH^ SCATTERING CATA 

The cross section for the formation of was measured at 
different angles of scatter by varying the detector distance Z 
while holding the magnet current I fixed. As pointed out in 
Section B of Chapter II the cross section is expected to be a pro- 
nounced peak at Hence the detector was swept through the 

range of scattering angles from 43^ to 49^. This was repeated for 
various target gas pressures for each value of the energy of the 
incident proton beam. The energies investigated were 70, 85, 100, 
150 and 200 electron volts. The scattering data recorded at each 
of these energies is shown in Figs. 26, 28, 30, 32 and 33 
respectively. A prominent peak exists in the 70, 85, and 100 eV 
data. The 150 and the 200 eV data shows no apparent peak in the 
cross section. 

The data was analyzed in the following manner. From the cross 
section measurements made at each energy an average value of CT ( 0 ) 
was computed and sketched in as a solid line. The assumed back- 
ground was then drawn as a dotted line. Subtracting out the back- 
ground from the 70, 85, aind 100 eV data, the points connected by 
the solid curve in Figs. 27, 29, and 3I were obtained. To find 
out what the actual peak in the scattering data looked like before 
it was detected by our detector of approximately 1^ resolution, 
the effect of the angular width of the detector was determined and 
the actual curve unfolded. The curves dotted on Figs. 27, 29, and 
31 represent the unfolded original peak in the scattering data. 
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-22 



Figure 26. Cross section for 

h "^ + CH ^ ^2 ^^3 

at 70 ev. (background 
included) 
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CT(0) X 10 ca' 
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-22 




These peaks are then the actual cross sections for the rearrange- 
ment reaction: 

+ CH, + CH^ 

- 4-23 

as the energy of the H varies from JO to 100 eV . All the peaks 
exhibit the same characteristics: a sharp leading edge at the lower 
angles, a flat top and a long decreasing tail extending out to 
about 51^- The peaks dependence on energy is exhibited in the 
following table: 











OJ 


CD 










P> 


p> 










w 


w 
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Mh cj 






r: 


H 


w U 


0 u 


0 


(t3 


4-> 


0 


0 






0) (/) 




'0 


•V 0) 


TJ (D 


d ^ 


Ct w 


•H 


c 


c 


'H 4-> 


4-> (T3 






p ^ 


D CD 


0 (D 


•H 0 


(P 0 




C3 


1 — 1 


W iH 


c cu 


(H 0 




(1) 


(TJ CL 


+-> CL 




u\ 0 


nJ 


(y CL 


CD a 


c a 


c3 ^ 


c 


Q) 


U 


LI 0 


0 0 


S 0 


< 


Cu 


< 


< u 


u u 



70 


eV 


-22 2 

5.95x10 cm 


45 . 15 ° 


1 . 7 ° 


10.2 


1 — " 

20.6 


45 . 87 ° 


00 


eV 


4 . 26 x 10 ^^cm^ 


45 . 1 ° 


1 . 5 ° 


6 • 56 


11. 56 


45 . 6 ° 


100 


eV 


-22 2 
3.0x10 cm 


46.25° 


1 . 5 ° 


4.62 


7.6 


46.7° 


S150 


eV 


-22 2 

^1,0 xlO cm 


Not 


Not 


Not 


Not 


Not 








Observed 


Obs er ved 


Obs er ved 


Obser ved 


Observed 



Table II, Summary of Experimental Data 

The peak does exhibit the expected physical properties. As the 
energy increases the magnitude of the peak decreases, the width 
narrows and the angle at which the peak occurs moves towards 46.9^. 
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pajB aqq, s^ubsb xdsj (d) ^^\xnz) ‘do:). :iapun lea aqq, 

XTgUl 
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sq^uasaxdax (^) aAxno *17^ ‘^xj ux uoxq.aas ssoxd xisDxq.a xoaqq. aqq, oq, 
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The energy dependence of the experimental cross section is compared 
to the theoretical cross section in Fig. 34- Curve (a) represents 
o(0)max> magnitude of the flat top. Curve (b) represents 

the area under the flat top. Curve (c) represents the total area 
under the curve. 
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Vo CONCLUSION 



The close agreement between the experimental data and the 
predictions of the ion molecule rearrangement theory proposed by 
Bates, Cook and Smith is evident in the following observations* 

(1) The sharp peak in the cross section predicted by Bates et.al., 

is observed in the 70? 85? and 100 eV data» The fact that the 

peak is not observed on the 150 and 200 eV data is not sur = 

prising for extrapolation of the experimental data to 150 eV shows 

"22 2 

that the magnitude of the cross section should be 1-0 x 10 cm . 
The 150 eV experimental data on Fig, 32 shows that a peak of this 
magnitude or smaller cannot be observed, 

(2) The theory predicted the peak to occur at The peaks 

observed were located between 45 ol^ and 46 , 25*^5 the latter cor- 
responding to the higher energy. As the accuracy of the approxi- 
mations made in the theory increases with higher energy, one 
expects the location of the peak in the experimental data to ap- 
proach 46 , 9 ^ as the energy is increased. This trend is observed 
exper imen tally , 

( 3 ) As is evident from Fig, 34 ? the experimental cross section has 
an energy dependence similar to that predicted by the theory, 

( 4 ) The cross section predicted by Bates et.al,, is an upper limit. 
Hence it is expected that the experimental cross section might be 
smaller than the theoretical cross section. The experimental cross 
section shown in Figc 34 is in fact about a factor of 3 O smaller 
than the theoretical upper limit. 



87 



(5) As some of the approximations made in the theory may no longer 
be valid at low energies one expects the peak to broaden as the 
energy deer eases » This trend is indeed observed in the 100, 85 , 
and 70 eV data shown in Figs, 3 I j 29 j and 27 respectively. 

It is felt that the facts mentioned above constitute a satis- 
factory verification of the ion molecule rearrangement theory of 
Bates, Cook and Smith as applied to the formation of in the 
r eact ion 

+ CH, hI + CH^ . 

- 4-23 
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APPENDIX I 



MULTIPLIER GAIN MEASUREMENT 

The following sketch shows hovj^ in principle, the gain measure- 
ment can be made. 




By recording the incident current I^ and the detector output I^ 

the gain can be computed from: 

I 



G = 



D 



To measure I^ the Dynode Strip was used as a Faraday Cup. Two 



problems then arose: 



8 , 



(1) The Dynode Strip, because of its large resistance (10 0) 
retains a large residual negative charge. Thus in using the Dynode 
Strip as a Faraday Cup, I^ had to be greater than 10 amps so 

as to insure that this residual negative charge did not adversely 
affect the measurement. 

-6 

(2) The Multiplier saturates when I^ ^ 5 x 10 amps. Thus we must 

have (5 X 10 )G. 

To overcome these two difficulties the gain measurement was 
performed in the following sequence of steps. 

Step 1. A beam of about 10 amps was focused on the Beam 
Collector . 
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step 2. The Beam Collector was swung away from the back of the 



Step 3 



Step 4 



Step 5 



Step 6 



Scattering Cell. This allowed the ion beam to fall 

directly on the multiplier. With F.S.I. and D.S.O. 

grounded, (see Fig. 17) amd a Keithley 6l0 Micro-Micro- 

ammeter connected to the D.S.I. this incident current 

* -10 

was measured to be I ^ 10 amps). 

. The Multiplier was now restored to the configuration shown 
in Fig. 17- The Dynode Voltage was reduced to approxi- 
mately 1,000 volts (this insured that the anode current to 
be measured in Step 4 was less than 10 ^ amps the satu- 
ration current of the Multiplier). 

f 

. The Multiplier was turned on and a detector current 

-7 

('^r; 10 amps) was measured at the anode. Thus, with the 
Multiplier settings as selected in Step 3 the gain was 




The incident bdam was now reduced by adjusting the 
focusing electrodes (see Fig. 5) a significant amount to 

f 

cause to decrease at least two orders of magnitude. 

IT — Q 

This new detector current (^ lO amps) was then re- 

D ^ ^ ^ n 

corded. Thus the incident beam must have been I. - —t — . 

1 G 

The Dynode Voltage was turned up to 1,500 volts. The 
Multiplier now had the configuration desired for the 
experimental measurements to be made later on. The de- 
tector signal now recorded allowed the gain of the 
Multiplier in this configuration to be computed from: 
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The above procedure was repeated at various values of detector 
distance Z and focusing magnet current I selected so as to give a 
measure of G versus the axial magnetic field B. The results of 
such measurements is shown in Fig. l8. 
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APPENDIX II 



The two computer programs SOLANG and CROSEC are listed in the 
following chapter . 

SOLANG, is the program which integrates the equation of the 



for predetermined values of the Focusing Magnet Current I, and 
the Incident Proton Energy E over a range of scattering angles Q. 
The integration is preformed using the D.H.P.G.C. Subroutine of 
the IBM System 36 O -67 Scientific Subroutine Package in the MAIN 
portion of the program, and in the two Subroutines FCT and OUTP . 
The resulting distances Z^, where the trajectory crosses the 
magnetic axis and the corresponding angles (see Fig. 21) are 
stored and later tabulated. 

The program now continues on to the Subroutine EXDMEGA and the 
succeeding subroutines, to compute the solid angle dQ( 0 ) for each 
of the trajectories. The solid angle is computed using 
equation (I 8 ) of Chapter IV, which is: 



A2^ has been determined by Bush, (Ref. 5) and is given by the 
following relations; (see Fig. (19) for explanation of symbols 
used) 

I for 32 . 2 ° >0 2: 26 . 6 ° 
o 



motion of the scattered ion (equations (^ and (9) of Chapter IV) 




AZ = (^■^■a) cotan 0 - 2b 

o ' ' o 
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II for 



26.6° ^02 17.3° 

o 

AZ = (h+a) cotan 0 - b 

o ' ^ o 



III for 17.3° s 0 s 0° 

o 



AZ = 2 h cotan 0 
o o 



The relation 



d0 

hz 



is determined by fitting the results of Z 



o ; E,I 



versus 0 that were obtained in the trajectory integrations by a 

^Zo 

polynomial Z = f( 0 ) . We then evaluate , ■ at the various 0 

O Q 0 

values. In this manner the solid angle dQ is computed for various 
values of 0. 

In order that the results may be tabulated for various Z values 
it is necessary to fit the results of 0 versus that were ob- 
tained in the trajectory integrations by a polynomial 0 = f(Z). 
Hence, on choosing a z, the corresponding 0 can be computed and 
in the manner described above, the solid angle can be determined. 

CROSfiC is the program which computes the cross section from 
the experimental measurements of the Detector Current and the 

Incident Proton Current I. which were recorded at various detector 

1 



distances Z for a particular proton energy methane pressure P 



sc 



and a particular Focusing Magnet Current B 

Because the data is recorded at various detector distances Z, 

the coefficients of the polynomial 0 = f(Z) that were determined in 

the program SOLANG are read into CROSEC on data cards and so allow 

the 0 corresponding to each Z to be determined. Rather than read 

into this program the solid angle at each point Z, the coefficients 

of the polynomial 0 ^ = that were determined in SOLANG are read 

in on data cards. Hence using: dO( 0 ) = 1 . 6 l 2 TTsin 0 /:|-^ ) AZ and 

o 
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the relations for AZ^ referred to in the discussion of the program 
SOLANG, we compute the solid angle at the Q corresponding to the Z 
where the data was recorded. Hence using equation (20) of Chapter 
IV: _ 

I .P G da(0) 

we compute the cross sections. 

Using the transformation from the Laboratory to the Center of 

. . + + 

Mass coordinates for the reaction H + CH, + CH^ we can de- 

- k -2 3 

termine „ and hence a(0)^ . . . The results are both 

C of M ' 'C of M 

tabulated and plotted. 
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